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Chemistry  Department 
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Salt  Lake  City,  UT  84112 


Abstract 


Anionic  states  of  Lia  were  studied  theoretically  using  a  full- 
valence  complete  active  space  (CAS)  self-consistent  field  (SCF) 
approach  to  explore  their  potential  energy  surfaces  followed  by 
quadratic  configuration  interaction  (QCI)  with  single,  double,  and 
approximate  triple  excitations  (SD(T))  calculations  to  more  accurately 
determine  relative  energies  and  electron  detachment  energies.  In 
addition  to  the  known  llg+  state,  one  triplet  (^A2')  and  two  quintet 
(5Ai"  and  5A2')  states  were  found  to  be  relatively  low-lying  and 
electronically  and  geometrically  stable  in  D3h  geometries.  All  of  these 
states  remain  electronically  stable  in  C2v  and  Dooh  geometries.  Hence, 
they  are  amenable  to  photoelectron  spectroscopy  detection. 
Preliminary  CAS  SCF  results  for  the  valence  isoelectronic  Naa*  are 
also  reported. 
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I.  Introduction 


Alkali  metal  trimers  may  be  considered  as  prototype  clusters 
involving  three  monovalent  s-orbital  radicals.  An  excellent  review  on 
the  neutral,  cationic,  and  anionic  alkali  metal  clusters  has  recently 
been  published  by  Bonacic-Koutecky  et  al.  [1].  Our  interests  in 
double-Rydberg  anions  [2]  as  well  as  in  chemical  bonds  between 
Rydberg  radicals  [3]  led  us  to  this  study  on  anionic  states  of  the  alkali 
metal  trimers. 

We  recently  studied  several  low  lying  doublet  and  quartet 
states  of  the  neutral  lithium  trimer  [4].  In  addition  to  the  well  known 
pseudo-rotating  ground  state  we  considered  excited  electronic  states 
which  were  recently  explored  using  the  two  photon  ionization 
technique  [5,6].  Our  results  suggested  a  new  assignment  of  the 
experimental  peaks. 

In  this  study  we  explore  anionic  states  of  the  lithium  and 
sodium  trimers.  Photodetachment  spectra  of  the  alkali  metal  trimers 
other  than  Lia*  were  reported  by  Bowen's  group  a  few  years  ago  [7], 
whereas  the  spectrum  of  Lia'  has  only  recently  been  recorded  [8]. 
The  lithium  clusters  are  difficult  for  experimental  studies  due  to  the 
high  temperatures  necessary  to  heat  the  metal  and  the  corrosive 
properties  of  lithium  [9]. 

Photoelectron  spectra  of  all  M3*  species  (M=alkali  metal  atom) 
are  dominated  by  two  distinct  and  well  separated  peaks.  These  peaks 
are  assigned  to  transitions  beginning  in  the  linear  (cyg2au2)  state 
of  M3*,  with  the  resulting  neutral  species  formed  either  in  the  ground 
^Zu'*'  (Og^duO  or  in  the  excited  2Zg+  (ag^Ou^)  state. 

For  the  ground  2Zu+  state,  the  vertical  detachment  energy 
(VDE)  is  1,309±0.005  and  1.158±0.010  eV,  for  Li3-  and  Na3', 
respectively  [7,8].  Upon  bending,  the  state  correlates  with  the 

lowest  2b 2  state  (a  local  minimum  on  the  pseudo-rotating  ground 
state  surface  of  the  neutral).  In  earlier  work,  we  indeed  verified  that 
linear  Li3  in  the  2Zu'*’  state  is  unstable^  with  respect  to  the  bending 
motion  [4].  Hence,  for  the  21^+  +  e  <-  lZg+  transition,  one  can  expect 
a  difference  between  the  VDE  and  adiabatic  electron  affinity  (EAa) 
equal  roughly  to  the  linearization  energy  of  the  ground  state  of  M3. 


For  the  second  +  e  <--  transition,  the  experimental 

VDE  is  2.003±0.005  and  1.813±0.010  eV,  for  Lis*  and  Nas*. 
respectively  [7,8].  Our  earlier  work  shows  that  the  ^Ig"^  state  is  stable 
with  respect  to  bending  motion  [4],  hence  the  discrepancy  between 
the  VDE  and  EAa  is  expected  to  be  small  for  this  transition. 

Past  theoretical  studies  on  the  anionic  alkali  metal  trimers 
concentrated  on  the  linear  structure.  Boustani  and  Koutecky  found 
the  state  of  Lia*  to  be  geometrically  stable  at  the  SCF  level  of 

theory,  which  was  consistent  with  the  electrostatic  repulsion 
argument  [10].  Electron  detachment  energies  from  the  linear 
structures  of  Li3-  and  Nas*  were  studied  by  Ortiz  using  electron 
propagator  theory  [11]  and  by  Bonacic-Koutecky  et  al.  using  the 
configuration  interaction  (Cl)  approach  [12].  These  theoretical  results 
match  the  experimental  data  with  an  accuracy  of  0.07  eV. 

In  the  present  work,  we  start  from  the  monomers  and  dimers 
of  the  alkali  metals  to  calibrate  the  accuracy  of  our  calculations.  The 
anions  of  M2  are  of  particular  interest.  The  photoelectron  spectrum  of 
Na2'  was  reported  in  Ref.  [7],  whereas  that  of  Li2‘  has  only  recently 
been  recorded  [13].  Three  observed  electronic  bands  are  assigned  to 
transitions  from  the  X2lu+  (Og^Ou)  bound  state  to  XiZg+  (Og^),  a^Zu^ 
(OgOu),  and  AlZu+  (<yg<Ju)  states  of  M2.  A  Franck-Condon  analysis  of 
the  lowest  electron  binding  energy  band  allows  one  to  determine 
spectroscopic  constants  of  7Li2*  [13]. 

Our  monomer  and  dimer  calibration  studies  are  followed  by  a 
thorough  exploration  of  the  singlet,  triplet,  and  quintet  states  of  Lis* 
for  D3h,  C2v,  and  Dooh  geometries.  To  the  best  of  our  knowledge,  the 

only  other  study  in  which  a  non-singlet  anionic  state  was  considered 
is  that  of  Boustani  and  Koutecky  [10].  They  concluded  that  a  linear 
stationary  point  of  the  lowest  triplet  state  is  not  a  minimum  at  the 
SCF  level  of  theory,  so  they  subsequently  focused  their  attention  on 
the  linear  singlet  state. 

In  our  study,  the  singlet,  triplet,  and  quintet  surfaces  are 
explored  at  the  full-valence  complete  active  space  (CAS)  self 
consistent  field  (SCF)  level  of  theory.  More  accurate  relative  energies 
and  electron  detachment  energies  are  determined  at  the  quadratic 
configuration  interaction  (QCI)  level  [14]. 
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Presentation  of  our  detailed  results  for  Li3'  is  followed  by  our 
preliminary  results  for  Na3'.  Guided  by  the  lithium  trimer  results, 
we  located  stationary  points  on  the  full-valence  CAS  SCF  potential 
energy  surfaces  of  the  neutral  and  anionic  Nas,  after  which  VDE’s  and 
EAa's  were  determined  at  the  same  level  of  theory. 

In  addition  to  the  lowest  energy  anionic  state,  we  have 

found  that  a  triangular  ^A2’  (ai'^e’^)  state  is  electronically  and 
geometrically  stable  for  both  alkali  metal  trimer  anions.  The 
electronic  stability  is  preserved  along  a  path  from  the  Dsh  minimum 
through  C2v  structures  to  a  linear  local  minimum  (producing  a 
state). 

Another  two  previously  unexplored  triangular  anionic  states  of 
5Ai"  (ai'e'2a2")  and  5A2'  (ai'e’^ai')  symmetry  were  also  found  to  be 
geometrically  stable.  They  are  also  electronically  stable  with  respect 
the  lowest  neutral  quartet  ‘^At*  (ai'e'2)  state.  Hence,  they  are  genuine 
bound  states  in  the  nonrelativistic  approximation  since  they  can  not 
decay  to  lower  lying  (e.g.,  Xl£g'*’)  anion  states  due  to  spin  and  orbital 

occupancy  restrictions.  The  electronic  stability  of  these  quintet  states 
is  preserved  along  a  C2v  path  which  connects  the  Dsh  and  Dc»h 

geometry.  At  linear  geometry,  these  two  quintet  states  correlate  to  a 
Sflg  (ogloulogljtu)  state.  The  5a  i"  state  preserves  a  positive 
curvature  along  the  antisymmetric  stretching  mode  in  C2  v 
geometries,  but  the  5  a  2'  state  (^82  in  C2v)  is  coupled  through  second- 
order  Jahn-Teller  interaction  with  another  quintet  5a  1  state  to  yield 
negative  curvature  along  the  b2  asymmetric  stretch  mode. 

In  Section  II  we  discuss  computational  aspects  of  this  study; 
Section  111  covers  the  anionic  monomers  and  dimers,  and  in  Section 
IV  our  results  for  the  singlet,  triplet,  and  quintet  states  of  Lis*  and 
Nas*  are  presented.  In  Section  V  we  summarize. 

II.  Computational  Aspects 

The  one-electron  atomic  basis  sets  for  the  lithium  atom  are 
detailed  in  Table  I.  The  smaller  one,  labeled  D,  consists  of  Dunning's 
double-zeta  (9s5p/3s2p)  one-electron  basis  set  [15]  supplemented 
with  diffuse  s  and  p  functions  with  the  same  exponent  0.0074  [16] 


and  one  d  function  with  the  exponent  0.2  [17].  The  larger  basis  set, 
labeled  T,  consists  of  Dunning's  triple-zeta  (I2s6p/4s3p)  set 
supplemented  with  the  diffuse  s  and  p  functions  [16]  and  two  d 
functions  with  exponents  0.3  and  0.15.  For  the  sodium  atom  we  used 
the  McLean  and  Chandler  (12s9p/6s5p)  set  [18]  supplemented  with 
diffuse  s  and  p  functions  with  the  same  exponent  0.0076  [19]  and 
one  d  polarization  function  with  the  exponent  0.175  [20]. 

The  ground  and  excited  potential  energy  surfaces  for  the 
neutral  and  anionic  M3  were  studied  at  the  CAS  SCF  level  of  theory, 
with  the  D  basis  set  for  Lis.  The  molecular  orbitals  derived  from 
atomic  core  orbitals  (Is  for  Li,  Is,  2s,  and  2p  for  Na)  were  doubly 
occupied  in  every  configuration  state  function  (CSF).  The  resulting 
neglect  of  the  core-core  and  core-valence  correlation  effects  yields 
results  we  label  frozen  core  approximation  in  subsequent  discussion. 
These  correlation  effects  are  negligible  for  lithium  and  small  for  the 
sodium  clusters  due  to  low  polarizability  of  the  M+  cores  [21].  The 
remaining  three  (neutral)  or  four  (anion)  electrons  were  distributed 
in  all  possible  ways  among  the  molecular  orbitals  derived  from  the 
atomic  valence  ns  and  np  orbitals.  This  choice  of  the  active  space  led 
to  572  and  220  CSF's  for  the  doublet  and  quartet  states  of  M3.  For  the 
anion,  the  number  of  CSF's  is  1716,  2145,  and  495  for  singlet,  triplet, 
and  quintet  states,  respectively.  The  CAS  SCF  calculation  for  the 
dimer  and  the  atomic  species  were  performed  with  the  analogous 
active  spaces  (i.e.,  frozen  core  approximation  plus  all  valence  orbitals 
in  the  active  space). 

The  CAS  SCF  calculations  were  performed  with  the  Utah  MESS- 
KIT  modular  electronic  structure  codes  [22]  which  generate  the 
analytical  second  geometrical  derivatives.  Geometrical  minima  and 
transition  states  on  the  three-dimensional  potential  energy  surfaces 
have  been  determined  using  our  automated  surface  walking 
algorithm  [23].  In  cases  where  a  small  energy  gap  between  electronic 
states  of  the  same  symmetry  occurs,  the  CAS  SCF  wavefunction 
optimization  for  the  higher  state  was  performed  in  the  state- 
averaged  (SA)  fashion  [24],  in  which  case  only  first  analytical 
derivatives  were  available  [25]. 


The  full-valence  CAS  SCF  method  is  sufficient  to  deliver 
detailed  information  about  the  potential  energy  surfaces.  In 
particular  it  can  deal  with  multiconfigurational  problems.  For  ^62 
anionic  state,  however,  the  nature  of  a  stationary  point  (minimum) 
determined  at  this  CAS  SCF  level  was  counterintuitive.  In  particular, 
the  curvature  along  the  asymmetric  stretching  mode  was 
questionable.  In  this  dubious  case,  we  also  examined  the  results 
obtained  using  CAS  active  spaces  smaller  and  larger  than  the  full- 
valence  space.  These  numerical  tests  strongly  suggest  that  the  ^62 
stationary  point  is  a  transition  state  rather  than  a  minimum. 

Due  to  the  small  number  of  valence  electrons  and  the  full- 
valence  active  space,  our  CAS  SCF  electron  detachment  energies  are 
quite  accurate.  For  instance,  the  VDE’s  for  Nas'  are  off  by  at  most 
0.08  eV,  for  the  two  experimentally  known  transitions.  Even  more 
accurate  detachment  energies  were  determined  for  Li3'  at  the 
quadratic  configuration  interaction  level  with  single,  double,  and 
approximate  triple  excitations  (QCISD(T))  [14].  This  method  takes  into 
account  dynamical  correlation  effects  and,  due  to  its  inherent  size- 
extensivity,  it  is  well  suited  to  deliver  reliable  electron  detachment 
energies.  Moreover,  in  our  QCI  calculations  we  removed  the  frozen 
core  approximation  so  all  core  electrons  were  correlated.  The 
QCISD(T)  results  were  obtained  with  the  Gaussian  92  suite  of  codes 
[26]  using  both  the  D  and  T  basis  sets.  The  difference  in  the  electron 
detachment  energies  did  not  exceed  0.04  eV,  thus,  to  conserve  space, 
only  the  T  basis  set  results  are  presented  in  later  discussion.  We 
expect  that  the  reported  detachment  energies  may  be  still  in  error 
by  ca.  0.02  eV  primarily  due  to:  (i)  inaccuracy  of  the  geometrical 
parameters,  (ii)  neglect  of  f-symmetry  basis  functions. 

Finally,  the  oscillator  strength  for  the  transition 

was  calculated  at  the  first-order  configuration  interaction  (FOCI) 
level  from  the  full-valence  CAS  reference  space  using  the  GAMESS 
package  [27]  and  the  D  basis  set. 
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III.  Atoms  and  Dimers 

Our  CAS  SCF  and  QCI  separation  between  the  atomic  2  s  and  2p 
states  of  Li  is  1.84  eV.  It  compares  well  with  the  experimental  data 
of  1.85  eV  [28].  For  Na,  however,  the  CAS  SCF  2S-2p  separation  is  too 
low  by  0.13  eV  in  comparison  with  the  experimental  value  of  2.10 
eV  [28].  This  larger  discrepancy  for  Na  is  caused  primarily  by  core¬ 
valence  correlation  effects  which  are  neglected  in  our  frozen-core 
approximation  [12]. 

Our  CAS  SCF  electron  affinities  for  Li  and  Na  are  0.52  and  0.47 
eV,  respectively,  compared  with  the  experimental  data  of  0.62  and 
0.55  eV  [29].  A  full-CI  calculation  for  the  two  outermost  valence 
electrons,  still  with  the  frozen  core  approximation,  leads  to  results  of 
0.59  and  0.54  eV.  Finally,  the  QCI  result  for  Li  with  the  frozen  core 
approximation  released  is  0.61  eV. 

These  results  demonstrate  that  our  full-valence  CAS  SCF 
approach  and  the  D  orbital  basis  set  are  capable  of  delivering  reliable 
S-P  energy  separations  between  valence  electronic  states  (accuracy 
ca.  0.15  eV)  as  well  as  meaningful  electron  affinities  (low  by  ca.  0.1 
eV).  The  inaccuracy  of  the  CAS  SCF  VDE's  is  primarily  caused  by 
insufficient  correlation  between  valence  electrons.  Refinement  of  the 
valence  electronic  states  separations  requires  inclusion  of  core¬ 
valence  and  core-core  correlation  effects. 

Our  CAS  SCF  equilibrium  internuclear  distances,  harmonic 
vibrational  frequencies,  and  dissociation  energies  for  the  four  lowest 
electronic  states  of  Li2  and  Na2  are  reported  in  Table  11.  The 
agreement  with  other  theoretical  and  experimental  results  is 
satisfactory.  The  discrepancies  in  equilibrium  geometries  do  not 
exceed  0.1  au,  with  exception  of  the  weakly  bound  a^Iu'*’  state  where 
due  to  the  shallow  potential  energy  depth  (De  «  0.02  eV),  a  larger 
discrepancy  (0.45  au)  is  reported  for  Na2.  The  harmonic  vibrational 
frequencies  match  the  experimental/theoretical  results  to  within  10 
cm‘1  which  is  better  than  expected.  The  CAS  SCF  and  QCI  errors  in  Do 
values  do  not  exceed  0.12  eV  and  0.04  eV,  respectively.  For  the 
excited  state  the  QCI  result  was  unreliable,  probably  due  to  the 
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serious  spin  contamination  at  the  UHF  level.  Only  the  CAS  SCF  results 
are  reported  for  this  state. 

In  Table  II  we  also  present  CAS  SCF  spectroscopic  parameters 
for  Li2'  and  Na2'.  The  agreement  with  the  available  experimental 
results  for  Li2'  [13]  is  excellent  for  Re  and  ©c,  but  D©  is  off  by  0.12 
(CAS)  and  0.04  eV  (QCI). 

In  Table  III  we  present  VDE’s  and  EAa’s  relative  to  the  bound 
22u'*'  ground  state  of  the  anion.  For  the  lowest  VDE  of  Li2‘,  the 
agreement  between  our  theoretical  (QCI)  and  the  experimental  value 
is  within  0.01  eV.  The  difference  between  the  VDE  and  EAa  is  the 
smallest  for  the  Li2(^Iu‘*')  +  e  <--  Li2'  transition  due  to  the  similarity 
in  equilibrium  geometries  of  the  neutral  and  the  anion.  The  Li2(^nu) 
+  e  <--  Li2*  transition  was  not  (yet)  observed  experimentally  [13], 
probably  due  to  its  two-electron  (i.e.,  shake  up)  nature.  The  CAS  SCF 
results  systematically  underestimate  the  VDE  by  ca.  0.1  eV  in 
comparison  with  the  experimental  results. 

IV.  Trimers 

Anionic  sates  of  Ms*  have  been  explored  in  Dsh,  C2v,  and  Dooh 
geometries.  The  dominant  electronic  configurations  of  the  numerous 
low-energy  anionic  and  neutral  states  can  best  be  clarified  using  the 
molecular  orbital  correlation  diagram  for  nsnp  atoms  shown  in  Fig.  1. 

Several  electronic  states,  (e.g.,  ^A2  ,  ^A  i",  5A2')  are 

electronically  stable  for  a  wide  range  of  molecular  geometries  and 
are  discussed  in  Sec.  IV  A.  Many  other  states,  such  as  ^Zg+,  ^Zu',  ^Zg+, 
^rig,  ^Ai',  25A2,  are  electronically  stable  only  for  a  limited  range  of 
molecular  geometries.  Certain  of  these,  such  as  the  25A2  state  are 
promising  candidates  for  long  lived  metastable  states.  Others,  such  as 
the  55:g+  state,  cross  the  energy  surface  of  the  ground  electronic  state 
of  the  same  multiplicity  and  thus  have  reduced  lifetimes.  All  of  these 
latter  electronic  states  will  be  discussed  in  Sec.  IV  B. 

The  dominant  electronic  configurations,  equilibrium 
geometries,  harmonic  vibrational  frequencies,  and  relative  energies 
are  reported  in  Table  IV  for  electronically  stable  (globally  or  locally) 
states  of  Li3'.  Calculated  values  of  the  VDE's  and  EAa's  for  various 


anion  states  are  shown  in  Tables  V  and  VI  for  Lis*  and  Na3-, 
respectively. 

A.  Globally  Electronically  Stable  Sates 
1.  D3h  Geonrivifies  and  the  C2v  Neighborhood 

The  lowest  energy  electronic  configuration  in  Dsh  geometry, 
(2ai’22e’2)  produces  the  following  electronic  states:  13A2’,  l^E',  and 
1 1 A 1 All  of  these  states  are  electronically  stable  with  respect  to  the 
neutral  ^E'(2ai'22e').  In  addition,  two  quintet  Lb*  anionic  states  were 
found  to  be  electronically  stable  with  respect  to  the  lowest  quartet 
neutral  state  ^A2'(2ai'2e’2).  The  latter  are  the  ^Ai"  (2ai'2e’2la2")  and 
5A2'  (2ai'2e'23ai')  states. 

a.  3A2’  state. 

This  is  the  lowest  energy  anionic  state  in  D3h  geometries.  The 
CAS  SCF  equilibrium  M-M  separations  are  5.93  and  7.09  au  for  Lia* 
and  Naa*,  respectively.  The  Lij*  harmonic  vibrational  frequencies 
(187  (e')  and  235  (ai*)  cm'l)  indicate  that  this  is  a  true  minimum  on 
the  potential  energy  surface.  For  Li3+  and  Na3+,  our  CAS  SCF 
equilibrium  M-M  distances  are  5.76  and  6.67  au,  respectively, 
somewhat  shorter  than  the  anionic  equilibrium  geometries. 

For  Li3-,  the  VDE  corresponding  to  the  transition  ^E’  +  e  <--  ^A2' 
is  0.761  eV  at  the  QCl  level  and  0.679  eV  at  the  CAS  SCF  level. 
Interestingly,  the  neutral  M3  thus  produced  would  be  in  the 
neighborhood  of  the  ^E'  conical  intersection.  Hence  this  triplet  anionic 
state  offers  the  experimental  opportunity  to  prepare  M3  neutral 
clusters  with  initial  geometries  otherwise  unavailable. 

Since  the  ground  state  of  the  M3  neutral  is  subject  to  first- 
order  Jahn-Teller  (FOJT)  distortion,  the  value  of  EAa  differs  from  VDE 
by  at  least  the  FOJT  stabilization  energy  ->  2b  2).  In  fact,  the  QCl 
difference  is  0.134  eV,  whereas  our  stabilization  energy  is  only  0.065 
eV. 


The  QCl  VDE  for  the  ^A2’  (ai’e'^)  +  e  <--  3A2’  transition  is  1.755 
eV  for  Li3‘.  The  resulting  neutral  '*A2'  is  the  lowest  energy  quartet 
for  D3h  geometries.  The  difference  between  this  VDE  and  EAa  is 
predicted  to  be  small  since  the  ^A2'  neutral  is  also  geometrically 
stable  at  D3h  geometries  with  an  equilibrium  M-M  distance  of  6.01 
au. 

It  is  worth  noting  that  our  electron  detachment  energies  from 
the  3 A  2'  state  discussed  above  (0.761  and  1.755  eV)  differ  from  the 
experimental  peaks  at  1.309  and  2.(X)3  eV  [8]  which  derive  from  the 
linear  singlet  form  of  Li3-.  This  makes  detection  of  the  3a 2’  anion 
plausible,  providing  a  significant  concentration  of  the  triplet  isomer 
can  be  produced  in  the  anionic  beam. 

For  today,  there  is  no  experimental  evidence  for  the  3 a  2’  state 
of  the  alkali  metal  trimers.  However,  an  unidentified  excited 
electronic  state  of  €03-  was  observed  by  Leopold  et  al.  in 
photoelectron  spectroscopy  experiments  using  2.540  and  2.707  eV 
photons  [30].  The  strong  peak  with  VDE  of  2.35-2.55  eV,  assigned  to 
a  linear  structure,  was  accompanied  by  a  feature  which  is  more  than 
two  orders  of  magnitude  weaker  and  possesses  a  VDE  of  ca,  1.5  eV.  It 
was  verified  that  the  weak  feature  is  due  neither  to  a  copper  oxide 
(e.g.,  Cu204')  nor  to  a  copper  cluster  containing  any  other  foreign 
species  [30].  Thus  the  authors  of  Ref.  [30]  assigned  the  weak  band  to 
an  electronically  excited  state  of  Cu3'.  Of  course,  the  peaks  observed 
for  CU3-  have  much  higher  VDE's  than  those  for  the  alkali  metal 
trimers,  consistent  with  the  difference  in  the  bulk  work  functions. 
Hence,  it  may  well  be  that  the  3  a  2'  state  of  Cu3*  was  observed  by 
Leopold  et  al.  [30]. 

b)  *E'  state 

The  state  l^E'  is  subject  to  FOJT  distortion  which  lowers  the 
molecular  symmetry  to  C2v  and  produces  two  distinct  electronic 
states 


l^Ai:  ci3ai22b22  -  C23ai24ai2 


11B2:  3ai24ai2b2 


with  the  HAiand  IIB2  states  favoring  obtuse  and  acute  triangle, 
respectively.  The  coefficients  ci  and  C2  are  equal  in  D3h  geometries, 
but  not  at  the  distorted  structures. 

The  UAi  state  correlates  with  the  lowest  energy  linear 
isomer  of  M3-  and  no  barrier  was  observed  along  this  rearrangement 
path.  As  M3-  approaches  the  linear  structure,  the  role  of  the  3ai24ai2 
configuration  decreases  (i.e.,  the  absolute  value  of  the  C2  coefficient 
uniformly  decreases). 

For  the  IIB2  offspring  of  the  UE’  state  the  full  valence  CAS  SCF 
geometry  optimization  led  to  a  stationary  point  with  R=6.193  au  and 
an  angle  of  55.19°.  The  harmonic  frequencies  of  all  normal  modes  are 
real  (see  the  seventh  row  in  Table  VII)  which  initially  led  us  to 
conclude  that  another  electronically  stable  minimum  for  M3-  had 
been  found. 

However,  when  we  performed  additional  CAS  SCF  calculations 
with  active  spaces  smaller  and  larger  than  the  full-valence  space,  the 
harmonic  frequencies  reported  in  Table  VII  were  obtained.  The 

frequencies  for  the  symmetric  stretching  and  bending  modes  are 
stable  within  35  cm-l.  The  asymmetric  stretching  mode  frequency  is 
imaginary  for  small  CAS  spaces,  real  for  the  full-valence,  and  again 
imaginary  for  the  largest  CAS  space  considered.  Clearly,  the  full- 

valence  CAS  SCF  prediction  as  to  the  nature  of  the  IB2  stationary 
point  seems  to  be  misleading.  This  is  an  annoying  artifact  of  non-full 
CI  calculations. 

The  behavior  of  the  b2  mode  frequency  as  a  function  of 

calculation  quality  must  be  related  to  symmetry -breaking  [31]. 
Previous  numerical  observations  indicated  that  this  may  happen  in 
single-configurational  calculations  in  which  some  correlation  effects 
could  be  included  in  the  wavefunction  at  the  expense  of 
incorporating  extra  configurations  of  the  wrong  symmetry.  A  remedy 
against  this  artifact  was  to  include  in  the  configuration  space  all 
important  configurations  of  the  correct  symmetry.  Here,  however. 


we  report  anomalies  at  the  full-valence  CAS  SCF  level  which  is 

believed  to  be  satisfactory  in  chemical  applications,  though 
computationally  demanding.  On  the  other  hand,  the  b2  mode 
frequencies  are  consistently  imaginary  in  small  valence  spaces.  The 
reported  results  indicate  that  the  anionic  l  B  2  stationary  point  is 
actually  a  transition  state. 

The  origin  of  the  b2  frequency  instability  definitely  requires 
further  investigations.  Our  preliminary  results  indicate  that  the  full- 
valence  CAS  SCF  energy  is  very  flat  with  respect  to  the  ai  <— >  b2 
orbital  rotations.  This  may  affect  the  stability  of  the  solutions  of  the 
coupled-MCSCF  equations  for  the  first-order  geometrical  responses. 
The  flatness  of  the  full-valence  energy  is  probably  related  to  the 
very  small  number  of  valence  electrons  in  the  case  of  alkali  metal 
clusters. 

The  three  equivalent  IB2  transition  states  connect  three 

equivalent  linear  ground  state  structures  of  symmetry,  with  the 

reaction  path  going  through  Cs  geometries.  The  graph  representing 
this  nonrigid  intramolecular  rearrangement  is  shown  in  Fig.  2.  The 
initial  distortion  from  the  linear  structure  involves  the  motion  along 
the  bending  (jtu)  and  antisymmetric  stretching  (Ou)  modes,  whereas 
that  from  the  C2v  structure  is  primarily  along  the  antisymmetric 

stretching  b2  mode.  The  energy  barrier  for  the  rearrangement  is  0.68 
eV  at  the  CAS  SCF  level  of  theory.  The  QCI  approach  is  inapplicable 
for  the  IB2  state,  similarly  as  for  the  state  of  the  dimer. 

c)  5  a  1"  state 

This  anionic  state  may  by  characterized  as  the  lowest  energy 

quartet  state  of  the  neutral  (i.e.,  '^A2’  (2ai'2e'2))  with  an  extra,  high- 
spin  coupled  electron  in  the  la2"  orbital  whose  electronic  density  is 
distributed  above  and  below  the  molecular  plane.  Hence,  it  does  not 
significantly  disturb  the  bonding  of  neutral  '^A2'.  for  which  the 
electron  density  is  accommodated  primarily  in  molecular  plane. 

Indeed,  the  equilibrium  R  for  this  5a  1"  state  (5.80  au)  is  even 
smaller  than  that  of  the  ^A2'  neutral,  which  is  6.01  au,  probably  due 
to  the  constructive  interaction  among  the  np  AO’s  in  the  la2"  MO.  The 
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similarity  between  the  neutral  and  the  anion  ^Ai"  manifests  also 
in  the  CAS  SCF  harmonic  frequencies  which  are  187  (e')  and  221  (ai') 
for  5a i"  and  189  (e')  and  218  (ai’)  for 

In  the  nonrelativistic  approximation,  the  ^Ai"  state  is  a  genuine 
bound  state.  The  VDE  to  the  lowest  neutral  state  is  0.56  eV 

which  is  even  lower  than  EAa  values  for  the  ^82  +  e  <--  ^kz 
transition.  Therefore,  the  5  a  1"  state  may  also  be  amenable  to 
photoelectron  spectroscopy  detection. 


d)  5a  2’  state 

This  state  may  be  characterized  as  the  neutral  ^A2'  with  an 
extra  high-spin  coupled  electron,  this  time  in  the  3ai'  orbital.  The 
extra  electron  density  is  in  the  molecular  plane,  so  larger  differences 
are  expected  between  the  neutral  and  the  anion,  than  those  reported 
above  for  the  5a  1"  ion. 

Indeed,  despite  the  bonding  properties  of  the  3a  1’  orbital,  the 
equilibrium  internuclear  distance  in  5a 2'  is  larger  by  0.3  au  than  in 
the  neutral  ^ki'  and  the  vibrational  frequencies  for  the  5A2’  state  are 
130  (e')  and  184  (ai'),  thus  smaller  than  those  for  the  ^A2'  neutral. 

In  the  nonrelativistic  approximation,  the  5 a 2'  anion  is  another 
bound  electronic  state.  The  transition  5a 2'  <--  5Ai"  has  a  vertical 
excitation  energy  of  994  cm*^  (QCI)  or  1178  cm'l  (FOCI)  and  the 
oscillator  strength  0.013  (FOCI).  This  prediction  should  help  in  future 
experimental  search  for  anionic  quintet  states. 

The  extra  electron  is  less  bound  in  this  5  a  2'  state  than  in  the 
5Ai"  state.  For  Lis'  the  value  of  the  5A2’  VDE  is  0.29  (CAS)  and  0.47 
(QCI).  For  Naa*,  the  VDE  value  of  0.12  eV  (CAS)  is  even  smaller. 

The  5a 2'  state  correlates  with  a  higher  in  energy  5ng  linear 
state,  as  the  5a  1"  state  does.  Both  5a 2'  and  5Ar'  states  may  be 
considered  as  molecular  analogs  of  the  core-excited  atomic  anionic 
states  studied  by  Bunge  and  collaborators  [32]. 


2.  Dooh  Geometries  and  the  C2v  Neighborhood 


a)  state 

This  is  the  lowest  energy  isomer  of  M3'.  The  Cl  wavefunction  is 
dominated  by  the  configuration  3og22ou^  and  the  2a u  orbital 
displays  an  antibonding  interaction  between  the  terminal  Li  atoms; 
hence  the  nature  of  the  doubly  occupied  2au  orbital  helps  preserve 
linear  structure.  The  same  effect  may  cause  the  equilibrium  M-M 
distance  to  be  larger  (by  0.24  au)  than  in  the  ground  state  linear 
neutral  22^+  (3ag22ou). 

The  latter  state,  however,  is  geometrically  unstable  with 
respect  to  bending  motion  since  the  effect  of  the  singly  occupied  2au 
orbital  is  apparently  overweighed  by  the  preference  of  the  doubly 
occupied  3og  orbital  for  Dsh  geometry.  The  2i„+  state  correlates  with 
the  2b  2  state  which  has  a  minimum  at  R=5.34  au  and  an  angle  of 
72.630. 

The  VDE  of  the  +  e  <--  lZg+  transition  is  1.322  eV  and  the 
EAa  is  smaller  by  0.317  eV,  whereas  the  linearization  energy  for  2B2 
neutral  is  0.289  eV.  Due  to  significant  difference  between  the  anionic 
and  neutral  equilibrium  structure,  the  photodetachment  peak  is 
expected  to  display  a  significant  Franck-Condon  width. 

Another  neutral  state  that  is  related  to  the  ^Zg'*'  anion  through 
a  single  electron  process  is  the  ^Zg'*'  state,  with  the  dominant 
configuration  3og2au^.  This  electronic  configuration  suggests  a  linear 
equilibrium  structure  and  a  bond  length  even  longer  than  in  the  iZg-*- 
anion.  This  is  indeed  the  case  as  the  equilibrium  R  for  ^^g'*'  is  5.96  au. 

The  photodetachment  peak  related  to  the  22;g+  +  e  <—  lZg+ 
transition  has  a  VDE  of  1.992  eV.  Since  there  is  small  relaxation  in 
the  neutral  equilibrium  geometry,  the  EAa  of  1.995  eV  is  quite 
similar  to  VDE,  so  the  peak  is  expected  to  be  sharp. 

Our  theoretical  results  are  in  a  good  agreement  with  the 
experimental  data  which  predict  VDE’s  of  1.309  ±  0.005  eV  (wide 
peak)  and  2.003  ±  0.005  eV  (sharp  peak)  for  transitions  leading  to 
2Zu'*’  and  22:g-f  products,  respectively  [8]% 


b)  ^Zu'*’  state 


The  dominant  configuration  of  this  state  is  3og22au4ag  which 
correlates  with  the  D3h  ^A2’  state.  Since  the  orbital  energies  of  2au 
and  4ag  have  opposite  dependence  on  bending  angle,  the  surface  of 
the  state  is  expected  to  be  quite  flat  with  respect  to  bending.  In 
fact,  our  bending  frequencies  are  only  53  cm-l  at  both  the  CAS  SCF 
and  MP2  levels.  Since  they  are  real,  a  transition  state  is  anticipated 
when  a  Dooh  structure  is  approached  from  the  D3h  minimum.  We 

actually  performed  a  transition  state  search  and  found  one  at  an 
angle  of  142°  and  R=5.95  au.,  with  an  energy  higher  by  0.271  eV 
(CAS)  than  at  the  D3h  minimum,  and  an  imaginary  frequency  of  34 
cm-l.  This  transition  state  was,  however,  0.(X)2  eV  lower  than  the 
linear  35;^+  minimum.  Approaching  the  linear  structure,  another  C2v 
minimum  was  found  at  R=S.99  au  and  an  angle  of  169o  with  an 
energy  only  0.002  eV  below  the  transition  state  and  a  real  bending 
frequency  of  23  cm*l.  We  abandoned  looking  for  other  transition 
states  since  the  energy  changes  involved  are  below  the  accuracy  of 
our  CAS  SCF  approach.  The  most  reasonable  conclusion  to  draw  is 
that  the  surface  is  actually  very  flat  with  respect  to  the  bending 
motion.  Interestingly,  the  optimal  distance  R  remains  within  the 
range  of  5. 9-6.0  au  for  all  the  stationary  points  found. 

The  32y+  state  is  electronically  stable  for  linear  geometries,  and 
the  value  of  EAa  for  the  transition  leading  to  the  lowest  energy  22:^+ 

neutral  is  still  0.589  eV.  We  thus  conclude  that  the  lowest  triplet 
anionic  state  is  electronically  stable  for  D3h,  C2v>  ^d  Dooh  geometries. 

c)  ^rig  state 

The  dominant  configuration  of  this  state  is  3og2Gu4agl}Cu- 
Despite  simple  molecular  orbital  arguments,  see  Fig.  1,  and  the  fact 
that  this  state  correlates  with  the  two  lowest  D3h^Ai"  and  ^A2’ 
states,  the  ^rig  state  lies  above  the  ^Zg+  (3ag2ou4ag3ou)  state  by 
0.117  eV.  The  lower  5ig+  state,  which  is  electronically  stable  only  in 
the  neighborhood  of  the  linear  structure,  will  be  discussed  in  Sec. 
IVB. 

The  CAS  SCF  frequencies  for  the  ^rig  state  indicate  that  the 
Renner-Teller  splitting  for  the  bending  modes  is  less  than  0.5  cm-l. 
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The  fact  that  both  such  frequencies  are  real  indicates  that  the  ^Ai" 
(5A2)  and  ^A2'(^B2)  states  correlate  to  the  ^flg  state  through 
transition  states.  Our  searches  for  these  transition  states  started  at 
D3h  geometries  and  moved  uphill  along  the  bending  mode.  Stationary 
points  at  R=6.02  au,  angle  126.8°,  vbcnd=  57i  cm"^  and  R=6.07  au, 
angle  131.4°,  vbend=  47i  cm'^  were  found  for  the  ^A2  and  ^82  states, 
respectively.  The  energies  of  these  transition  states  are  0.02  eV 
below  (5A2)  and  0.01  eV  above  (^82)  the  energy  of  the  ^flg  minimum. 
Clearly,  the  surfaces  become  very  flat  from  the  angle  of  ca.  120°;  for 
this  reason  we  did  not  pursue  further  study  for  C2v  geometries  closer 
to  the  linear  structure  since  the  accuracy  of  our  CAS  SCF  data  would 
be  not  sufficient.  The  flatness  of  the  ^A2  and  ^82  surfaces  is  probably 
related  to  the  opposite  angular  dependence  of  the  2b2  and  4ai  orbital 
energies,  see  Fig.  1. 

The  fact  that  the  5lg+  state  is  lower  than  the  ^rig  state  at  linear 
geometries  indicates  that,  in  a  C2v  geometry,  there  is  a  crossing 
between  the  ^Ai  (5Xg+)  and  the  ^82  component  of  the  ^flg  state.  In 
this  region,  the  ^82  state  becomes  geometrically  unstable  with 
respect  to  the  b2  vibrational  mode  due  to  second -order  Jahn-Teller 
coupling  [33].  In  addition,  "hopping"  to  the  ^Ai  (^Ig+)  surface  can 
take  place.  80th  factors  will  decrease  the  lifetime  of  the  5a  2'  state,  in 
addition  to  the  radiative  decay  channel. 

The  electronic  stability  of  the  5ng  state  has  to  be  considered 
with  respect  to  two  low-energy  neutral  quartet  states, 
(3ag2au4ag)  and  4ng  (3ag2oulJCu)-  The  equilibrium  geometries  of  the 
former  (R=7.47  au)  and  the  latter  (R=5.41  au)  differ  significantly 
from  the  equilibrium  R  of  ^iig,  which  is  6.08  au.  Hence,  a  major 
difference  is  expected  between  the  corresponding  VDE  and  EAa 
energies:  VDE=0.656  and  0.783  eV,  EAa  =  0.603  and  0.640  eV  for  the 
'♦Zu'*’  and  ^rig  transitions,  respectively. 

We  conclude  that  the  5a  1"  and  5A2'  quintet  anionic  states 
are  electronically  stable  in  Dsh,  C2v.  and  Dooh  geometries. 

8.  Locally  Electronically  Stable  Species 


a)  l^Ai'  state 


This  state  results  from  the  same  2ai’22e’2  configuration  that 
produces  the  lE'  and  ^\2  states.  Even  though  this  state  is  not  subject 
to  FOJT  distortion,  we  write  its  approximate  Cl  wavefunction 
ci3ai22b22  +  C23ai24ai2  in  the  C2v  orbital  symmetry  notation  to  make 
the  difference  between  UE'  and  l*Ai‘  transparent. 

The  geometry  of  this  state  was  optimized  within  Dsh  symmetry 
and  the  geometrical  hessian  analysis  indicates  that  this  is  a  minimum 
structure  indeed,  see  Table  IV.  The  anion  is  electronically  stable  with 
respect  to  the  lowest  ^E'  (2ar22e')  state  of  the  neutral  by  0.065  eV 
(CAS  level). 

Recall  that  the  ^E'  state  of  the  neutral  is  FOJT  unstable  and  we 
found  that  the  2b  2  minimum  of  Lia  is  lower  by  0.058  eV  (CAS)  than 
the  D3h  *Ai'  minimum  of  Lis*.  Clearly,  the  anionic  and  neutral 
surfaces  must  cross  in  C2v  symmetry  not  far  from  the  Dsh  manifold. 
The  electronic  stability  of  the  anion  is  thus  restricted  to  a  small  part 
of  the  potential  energy  surface  so  this  state  cannot  be  long  lived. 

b)  3Zg+  state 

This  is  a  state  with  two  important  configurations:  3ag2Gu^4ag 
and  3og22ou3ou-  Its  minimum  energy  geometry  was  found  at  R=6.04 
au  where  the  bending  mode  frequency  is  only  86  cm*l.  Hence,  C2v 
geometries  are  easily  accessible.  This  state  correlates  at  D3  h 
geometries  to  an  electronically  unstable  ^E'  state.  Therefore,  the 
electronic  stability  is  restricted  to  linear  geometries  and  their  close 
neighborhood;  at  linear  geometry,  the  state  is  electronically  stable 
with  respect  to  the  lowest  neutral  2z„+  state  by  0.340  eV  (EAa). 


c)  ^rig  state 

For  this  anionic  state,  the  .dominant  configuration  is 
3og22au37Cu.  A  geometry  minimum  was  found  at  R=5.53  au  with  a 
Renner-Teller  splitting  of  7  cm**.  The  bending  frequencies  are  low, 
hence  C2v  geometries  are  easily  accessible.  The  ^82  component  of  this 


^rig  state  correlates  with  the  same  D3h  ^E’  state  as  the  state 

does.  The  3^2  component  of  this  state  correlates  with  an 
electronically  unstable  state.  Therefore,  the  electronic  stability  of 
the  ^Ilg  state  is  restricted  to  linear  geometries  and  their  close 
neighborhood  where  it  is  stable  with  respect  to  the  lowest  neutral 

state  by  0.173  eV  (EAa). 

d)  52:g+  state 

This  is  the  lowest  energy  quintet  state  at  linear  geometries  and 
has  a  dominant  configuration  3og2ou4og3  Ou.  Due  to  the 
counterbalancing  energy  dependence  of  the  2au  and  4ag  orbitals,  the 
surface  is  expected  to  be  flat  with  respect  to  bending.  In  fact,  the 
bending  frequency  is  only  76  cm-l  and  Civ  geometries  are  easily 
accessible. 

In  Dah  symmetry,  this  state  correlates  with  the  ^E*  (2ai'2e'23e') 
state  which  is  electronically  unstable.  Hence,  the  electronic  stability 
of  the  ^Ig'*'  state  is  restricted  to  linear  geometries  and  their  close 
neighborhood,  where  a  crossing  with  the  (^Ai’)  state  takes  place, 
accompanied  by  second-order  Jahn-Teller  instability  and  possibility 
of  surface  "hopping". 

At  linear  geometries,  the  5ig+  state  is  electronically  stable  with 
respect  to  the  two  lowest  quartet  states  of  the  neutral:  VDE=0.726  eV 
and  1.073  eV,  EAa=0.720  eV  and  0.756  eV,  for  the  41^+  and  ^rig  state, 
respectively. 

e)  5e"  state 

The  5e"  state  is  a  candidate  for  a  long-lived,  metastable, 
pseudo-rotating  state.  The  electronic  stability  of  this  state  is 
discussed  below,  for  different  regions  of  the  potential  energy  surface. 

In  D3h  geometries,  the  (2ai '2e' 1  a2 " 3 a i ')  state  is 

electronically  unstable  with  respect  to  the  ground  quartet 
(2ai’2e'2)  state.  However,  since  the  2ai’2e' la2"3ai'  configuration 
differs  from  the  2ai'2e'2ke"  configuration  needed  to  produce  the  4a 2' 
neutral  (ke"  denotes  an  e"  symmetry  orbital  for  the  outgoing 


electron)  by  two  spinorbital  occupancies,  the  coupling  to  this 
detachment  channel  is  expected  to  be  weak.  Among  configurations 
which  contribute  to  ^E"  and  simultaneously  can  couple  to  the  ‘*A2 
neutral  in  a  one-electron  process,  the  largest  Cl  coefficient  (only 
0.12)  was  found  for  the  2ar2e'2le"  configuration.  The  neutral  state 
(2ar2e'la2").  to  which  the  anion  could  couple  in  a  one-electron 
process,  is  higher  in  energy  than  the  anion  by  0.23  eV  (CAS  level  at 
the  cation  equilibrium  geometry). 

The  5E"  state  is  FOJT  unstable,  producing  2^A2  and  l^Bi  states 
which  favor  obtuse  and  acute  structures,  respectively,  consistent 
with  the  bending  angular  dependence  of  the  4a i  and  2b2  orbitals.  A 
stationary  point  on  the  (3ai4ai  IbiSai)  surface,  with  R=6.44  au 
and  an  angle  of  43. 6^,  proved  to  be  a  transition  state,  with  a  negative 
curvature  along  the  asymmetric  stretching  mode,  see  Table  IV.  At 
this  geometry,  the  energy  of  the  5Bi  state  is  lower  than  the  energy  of 
the  lowest  quartet  ^B2  neutral  state  by  0.04  eV  (CAS)  and  0.22  eV 
(QCI).  Hence,  the  electronic  instability  of  the  anionic  quintet  state  is 
probably  restricted  to  Dsh  geometries  and  the  close  neighborhood. 

Geometry  optimization  for  the  2^A2  offspring  of  the  ^E"  state 
was  performed  in  the  state-averaged  CAS  SCF  procedure,  to  avoid 
"root-flipping"  [34].  Three  electronic  states  were  state-averaged  with 
equal  weighting  factors:  15A2(^Ai"),  1^B2(^A2’),  and  25A2  (^E"),  and 
the  Cl  mixing  was  significant  between  two  ^A2  3ai2b24ailbi  and 
3ai2b2lbi5ai  configurations. 

A  stationary  point  for  the  25 A 2  state  was  found  at  R=5.28  au, 
and  an  angle  of  111.5°.  We  expect  it  is  a  minimum,  but  the  finite 
difference  hessian  formed  in  the  SA  CAS  SCF  C2v  optimization  does 
not  allow  us  to  determine  the  b2  mode  frequency.  Interestingly,  at 
this  stationary  point,  the  energy  of  the  lowest  4b  2  neutral  state 
calculated  without  the  SA  constraint  is  higher  by  0.14  eV  than  the 
energy  of  the  SA  optimized  25  A  2  anion.  This  indicates  that  the  25  A  2 
state  is  electronically  stable  at  the  stationary  point. 

Determination  of  the  pseudo  rotation  barrier  on  the  5e"  surface 
requires  a  consistent  calculation  of  the  15Bi  and  25 A  2  energies. 
Hence,  we  recalculated  the  energy  of  15B  1  at  its  stationary  point 
using  the  SA  procedure  in’ olving  15A2  ,  1^B2  ,  and  15Bi  states.  The 


resulting  energy  was  only  0.07  eV  higher  than  the  single-state  CAS 
SCF  energy  and  the  pseudo  rotation  barrier  is  estimated  to  be  544 
cm-l. 

The  2^K2  state  correlates  with  the  linear  (3ag2aul tcu^) 
state.  This  i'  the  third  electronically  stable  anionic  quintet  state  at 
the  linear  geometry,  and  is  higher  in  energy  than  the  and  ^Flg 

states  by  0.341  and  0.225  eV  (QCI),  respectively.  A  stationary  point 
was  found  for  the  ^lu'  state  at  R=5.29  au,  with  a  bending  frequency 
of  only  56  cm-l.  The  anion  is  electronically  stable  with  respect  to  the 
two  lowest  quartet  states  of  the  neutral:  VDE=0.678  and  0.410  eV, 
EAa  =  0.379  and  0.415  eV,  for  and  '^Ilg,  respectively. 

V.  Summary 

The  alkali  metal  trimer  anions  possess  more  than  one 
electronically  and  geometrically  bound  state.  In  addition  to  the 
known  ^2g+  state,  one  triplet  3A2'  and  two  quintet  ^Ai",  5a2’  states 
are  predicted  to  be  electronically  stable  for  a  wide  variety  of 
geometries.  The  multitude  of  electronically  bound  states  is  probably 
related  to  the  small  energy  separation  between  valence  ns  and  np 
atomic  orbitals  in  alkali  metal  atoms. 

The  linear  l2g+  (3og22ou^)  state  is  the  lowest  energy  isomer  of 
Li3-.  Our  VDE  for  the  22^+  +  e  <-  ^Ig+  and  22g+  e  <-  i2g+ 
transitions  are  1.322  and  1.995  eV,  respectively,  in  good  agreement 
with  the  experimentally  observed  peaks  at  1.309  ±  0.005  and  2.003 
±  0.005  eV  [81.  The  barrier  for  nonrigid  intramolecular 
rearrangement  of  the  linear  isomer  was  found  to  be  0.68  eV  and  to 
pass  through  a  C2v  symmetry  ^82  transition  state. 

The  3  A  2'  (ai  '2e'2)  anion  is  the  lowest  energy  isomer  in  D3  h 

geometries.  This  triplet  state  remains  electronically  stable  in  C2v  and 
Dooh  geometries.  Our  VDE  for  the  2e'  +  e  <—  ^A2'  and  ^A2'  +  e  <— 

transitions  in  Li3*  are  0.761  and  1.755  eV,  respectively.  Hence,  the 
peaks  should  be  clearly  distinct  from  ,the  transitions  involving  the 
linear  lZg+  anion,  thereby  suggesting  further  experimental  study. 

After  this  work  was  completed,  we  were  pleased  to  find  that 
Bauschlicher  et  al.  [35]  also  invoked  the  ^A2'  anionic  state  to  explain 


the  weak  feature  in  the  photoelectron  spectrum  of  Cus*  [30]  and 
found  a  VDE  for  the  triplet  Dbh  isomer  of  CU3-,  calculated  at  the 
modified  coupled  pair  functional  level  of  theory,  in  good  agreement 
with  the  experimental  finding.  Interestingly,  the  weak  feature 
ascribed  to  the  3a  2'  anion  was  not  observed  in  the  later  work  from 
the  Lineberger  group  [36],  presumably  due  to  cooler  ion  source 
conditions. 

Two  quintet  ^Ai"  and  5 A  2’  anionic  states  were  also  identified 
that  may  be  thought  of  as  the  lowest  energy  quartet  ^A2  (ai'e'2) 
state  of  the  neutral  with  an  extra  high-spin  coupled  electron  attached 
to  a  a2"  or  ai'  orbital,  respectively.  These  quintet  states  are  genuine 
bound  states  in  the  nonrelativistic  approximation,  and  may  be 
considered  as  molecular  analogs  of  core-excited  atomic  anionic  states 
studied  earlier  [32].  The  vertical  excitation  energy  for  the  radiative 
5A2'  <"  ^Ai"  transition  is  ca.  1000  cm*I  and  the  oscillator  strength  is 
0.013. 

For  the  ^Ai"  and  5 A2’ states  of  Lia*,  the  vertical  detachment 
energies  to  the  ^^2  state  neutral  are  0.563  and  0.465  eV, 
respectively.  Both  quintet  states  remain  electronically  stable  with 
respect  to  the  lowest  neutral  quartet  state  along  a  bending  motion 
path  to  the  linear  structure,  where  they  correlate  with  a  ^rig  state. 
When  proceeding  to  the  linear  structure,  the  5a2'(5B2)  state 
intersects  with  the  52:g+(5Ai)  state,  so  geometrical  instability  must 
develop  along  the  b2  mode  and  surface  "hopping"  may  take  place. 

The  lithium  trimer  also  possesses  many  anionic  states  which 
are  electronically  stable  for  restricted  geometries  only.  The  most 
interesting  of  them  is  probably  the  pseudo-rotating  5e"  state.  It  is 
electronically  unstable  in  Dsh  geometries,  but  stable  at  the  stationary 
points  of  the  (transition  state)  and  5A2  (minimum)  surfaces. 

Moreover,  the  5e"  wavefunction  displays  properties  required  to  be 
classified  as  a  Feshbach  resonance  state  at  Dsh  geometries. 

An  extension  of  our  present  conclusions  to  other  alkali  metal 
trimers  would  be  important.  It  is  encouraging  that  our  preliminary 
CAS  SCF  results  for  Naa*  are  qualitatively  the  same  as  those  for  Liy. 
The  presence  of  alkali  metal  atoms  in  heteroatomic  radicals  may 
enhance  an  ability  to  form  more  than  one  electronically  bound 


anionic  state.  Indeed,  the  whole  family  of  lithium  substituted  double- 
Rydberg  anions  accommodates  an  excited  triplet  state  in  addition  to 
the  ground  fully-symmetric  singlet  state  [37]. 
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Table  I.  The  Gaussian  one-electron  basis  sets  used  for  the  lithium 
atom  [15-17]. 
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symmetry  exponent  contraction  coefficients 

D  basis  set 


s 

1469.0 

0.000766 

-0.000120 

220.5 

0.005892 

-0.000923 

50.26 

0.029671 

-0.004689 

.  14.24 

0.109180 

-0.017682 

4.581 

0.282789 

-0.048902 

1.580 

0.453123 

-0.096009 

0.564 

0.274774 

-0.136380 

0.07345 

0.009751 

0.575102 

0.02805 

-0.00318 

0.517661 

s 

0.02805 

1.0 

s 

0.0074 

1.0 

p 

3.2500 

0.008691 

0.6474 

0.047819 

0.1689 

0.210615 

0.05561 

0.528695 

0.02047 

0.381375 

p 

0.02047 

1.0 

p 

0.00740 

1.0 

d 

0.2 

1.0 

T  basis  set 

s 

6601.0 

0.000117 

-0.000018 

989.7 

0.000911 

-0.000142 

225.7 

0.004728 

-0.000741 

64.29 

0.019197 

-0.003020 

21.18 

0.063047 

-0.010123 

7.724 

0.163208 

-0.027094 

3.003 

0.314827 

-0.057359 
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1.212 

0.393936 

0.493 

0.196918 

0.09515 

0.009997 

0.04791 

-0.005402 

0.02220 

0.001704 

s 

0.04791 

1.0 

s 

0.02220 

1.0 

s 

0.00740 

1.0 

p 

5.928 

0.003668 

1.295 

0.020596 

0.3493 

0.083282 

0.1154 

0.279201 

0.04401 

0.514145 

0.01784 

0.278900 

p 

0.04401 

1.0 

p 

0.01784 

1.0 

p 

0.0074 

1.0 

d 

0.3 

1.0 

d 

0.15 

1.0 

-0.093895 

-0.121091 

0.276608 

0.549548 

0.277385 


Table  11.  Full  valence  CAS  SCF  spectroscopic  parameters  for  different 
electronic  states  of  M2  and  the  ground  state  of  M2’  (M=Li,Na).  R, 
in  au,  cDe  in  cm**.  Do  in  eV.  The  QCl  SD(T)  values  for  Do  are  in 
parenthesis. 


Species 

State 

Li2 

x>v 

Li2 

a3i:u+ 

Li2 

b3nu 

Li2 

AlXu-" 

Li2* 

X2Zu+ 

5 

5 


8 

0 


5.01 

4.96C 


6.01 

5.93c 


5.85 

5.85d 


©c 

Do 

343 

0.92  (1.02) 

35ia 

1.05a 

52 

0.02  (0.04) 

61*> 

0.03b 

335 

1.31  (1.46) 

342C 

1.42C 

256 

1.09 

261C 

1.14c 

230 

0.75  (0.83) 

232d 

0.87d 

Na2 

Xllg* 

Na2 

a3Zu+ 

Na2 

b3nu 

Na2 

Al2u+ 

Na2' 

X2Zu+ 

5.89 

5.82a 


10.29 

9.84e 


6. 

6. 


7,05 

6.88a 


6.85 


53 

59a 


20 

25f 


112 

117a 


98 


.66 

.72a 


0.02 

0.02® 


45 

1.00 

55f 

1.08C 

0.95 

l.OOa 


.55 

.44a 


a)Ref  [38],  ORef  [39],  c)Ref  [40],  b)Ref  [41],  d)Ref  [13],  c)Ref  [42] 
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Table  III.  CAS  SCF  and  QCl  SD(T)  vertical  electron  detachment 
energies  (VDE)  and  adiabatic  electron  affinities  (EAa)  in  eV  for  the 
state  of  X2'  (X=Li,  Na).  Experimental  results  in  parenthesis. 


Transition 

Method 

VDE 

EAa 

Li2(^Za+)-i-e  4-  Li2‘ 

CAS 

0.42 

0.34 

Li2(3Zu+)+e  <- Li2*  I  CAS 


Li2(^nu)+e  4—  Li2*  I  CAS 


1.37 

1.26 

1.497 

1.406 

1.90 

1.79 

1.948 

1.807 

Li2(^Iu''')+e  4-  Li2' 

CAS 

2.02 

2.02 

Na2(^2e+)+e  f-  Na2'  I  CAS 


0.43  (0.543)b 


Na2(^Zu'^)+c  4-  Na2* 

CAS  J 

1 1.13 

1.00 

Na2(^nu)+e  4-  Na2'  j 

CAS 

12.04 

1.98 

Na2(li:u'^)+e  4-  Na2-  I  CAS 


2.04 


2.03 


aRef.  [13] 
bRef.  [7] 
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Table  IV.  Dominant  electronic  configurations  followed  by  the  Cl 
coefficients,  geometries  (distances  in  au,  angles  in  degrees),  CAS  SCF 
harmonic  frequencies  (in  cm'l)  and  relative  CAS  SCF  and  QCI  SD(T) 
energies  (in  eV)  for  stationary  points  on  different  potential  energy 
surfaces  of  Lis*. 


State 

dominant 

configuration 

Geometry 

Frequencies 

AEcas 

AEga 

'V 

3og22ou-  (0.84) 

Dooh 

R=5.757 

106(jcu) 

184(ag) 

284(au) 

0.0 

0.0 

lB2a) 

3ai24ai2b2 

(0.87) 

C2v  R=6.193 
Z=55.193 

173(aibend) 

252(ai) 

658i(b2) 

0.68 

3A2' 

2ai’22e'2  (0.88) 

D3h  R=5.926 

187(e’) 

235(ar) 

0.46 

0.379 

3og22ou4og 

(0.84) 

Dooh  R=6.030 

53(jtu) 

153(crg) 

178((yu) 

0.74 

0.705 

SAi" 

2ai'2e'2la2" 

(0.96) 

D3h  R=5.802 

187(e’) 

221(ar) 

1.70 

1.571 

5A2' 

2ai’2e’23ai' 

(0.97) 

D3h  R=6.311 

130(e') 

184(a,’) 

1.82 

1.700 

^ng 

3og2ou4og  1 

Dooh  R=6.083 

51(7Cu) 

135(Og) 

204(Ou) 

1.98 

1.902 

lAr 

2ar22e'2  (0.72) 

D3h  R=5.928 

245(ar) 

301(e’) 

1.08 

3og2ou24og 
(0.63)  + 
3og22ou230u 
(0.60) 

Dooh  R~6.038 

86(7Cu) 

161(og) 

303(au) 

1.02 

0.954 

^ng 

3<Tg^2cru  1  Jtu 
(0.77) 

DoohR=5.529 

68(jtu)->3A2 

75(jtu)->^R2 

184(Og) 

293(au) 

1.20 

1.122 

3og2ou^®g3®u 

(0.99) 

Dooh  R=6.558 

76(jiu) 

132(Og) 

261(au) 

1.81 

1.785 

5Bi 

3ai4ailbi5ai 

(0.95) 

C2v  R=6.439 
Z=43.640 

155(aibend) 

340(ai) 

100i(b2) 

2.29 

2.122 

5A2b) 

3ai2b2lbi5ai 
(0.84)+ 
3ai2b24ai  Ibi 
(0.41) 

C2v  R=5.276 
Z=lll.531 

27(aibend) 

273(ai) 

?(b2) 

2.22 

2.054 

5lu* 

3  0g20u  1 
(0.91) 

Dooh  R=5.289 

56(ffu) 

I92(ag) 

314(Ou) 

2.21 

2.126 

2)CAS  SCF  frequencies  obtained  with  the  active  space  of  seven  ai,  two 
a2,  three  bi,  and  five  b2  orbitals, 

l’)SA  CAS  SCF  geometry  optimization;  relative  energies  based  on  the 
5Bi  results  and  the  SA  CAS  SCF  difference  between  and  ^A2 


Table  V.  Electron  vertical  detachment  energies  (VDE)  and  adiabatic 
electron  affinities  (EAa)  in  eV  for  different  anionic  states  of  Li3-, 
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Transition 

Method 

VDE 

EAa 

2Zu++e  ♦-  lZe+ 

CAS 

1.27 

1.02a 

22u++e  <- 

oa 

1.322b 

1.005a 

2Ze++  e 

CAS 

1.94 

1.93 

2Za++  e  <-  lZg+ 

oa 

1.992C 

1.995 

2E'+e  <-  3A2' 

CAS 

0.68 

0.56a 

2E'+e  4-  3A2' 

oa 

0.761 

0.627a 

'^A2  +e  4-  3A2' 

CAS 

1.63 

1.63 

'*A2'+e  4-  3A2’ 

oa 

1.755 

1.759 

‘*A2+e  4-  3Ai" 

CAS 

0.40 

0.40 

^A2+e  4—  ^Ai" 

oa 

0.563 

0.567 

4A2'+e  4-  5A2' 

CAS 

0.29 

0.27 

^A2'+e  4-  5A2' 

oa 

0.465 

0.438 

a) the  state  of  the  neutral  is  ^62  for  the  adiabatic  transition 

b) The  experimental  VDE  is  1.309±0.005  eV  [8] 

c) The  experimental  VDE  is  2.003±0.005  eV  [8] 


> 
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Table  VI.  The  CAS  SCF  electron  vertical  detachment  energies  (VDE) 
and  adiabatic  electron  affinities  (EAa)  in  eV  for  different  anionic 
states  of  Nas’ 


Transition 

VDE 

EAa 

2Zu++e  <- 

1.13a 

1.05b 

QIS&SIQI 

1.73c 

1.70 

2E'+e  <-  3A2' 

0.78 

0.62b 

'*A2'+e  <-  3A2’ 

1.52 

1.44 

^A2+e  4—  ^Ai” 

0.25 

0.22 

4A2'+e  «-  5A2' 

0.12 

0.11 

®)The  experimental  VDE  is  1.158±0.01  eV  [7] 

b) The  state  of  the  neutral  is  ^82  for  the  adiabatic  transition 

c) The  experimental  VDE  is  1.813±0.01  eV  [7] 
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Table  VII.  Harmonic  vibrational  frequencies  (cm-l)  for  ^62  state  of 
Lis*  obtained  with  different  CAS  spaces.  The  full-valence  CAS  SCF 
geometry  was  used  (R=6.193  au,  angle  =  55.190).  Orbitals  lai,  2ai, 
and  lb2  are  doubly  occupied  in  every  configuration  state  function. 


Active  Space 

Vbending 

Vsym.  stretch 

Vasym.  stretch 

two  ai  one  b2 

184 

265 

227i 

two  ai  one  bi 
one  b2 

173 

256 

227i 

three  ai  one  bi 
one  b2 

154 

230 

368i 

four  ai  one  bi 
two  b2 

179 

257 

242i 

four  ai  one  a2 
two  bi  two  b2 

178 

256 

484i 

four  ai  one  a2 
two  bi  three  b2 

178 

254 

328i 

five  ai  one  a2 
two  bi  four  b2 
(full-valence) 

175 

254 

427 

six  ai  one  a2 
two  bi  four  b2 

177 

254 

1033 

six  ai  two  a2 
three  bi  four  b2 

176 

254 

1571 

seven  ai  two  a2  173  252  65 8i 

three  bi  five  b2 _ 
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Captions  for  Figures. 

Figure  1.  Molecular  orbital  correlation  diagram  for  the  lithium  trimer 
as  a  function  of  the  bending  angle  6  which  is  equal  to  60°  and  180° 
for  the  D3h  and  Dooh  structure,  respectively.  The  SCF  orbital  energies 

were  obtained  for  the  closed-shell  cation  Li3+,  with  the  3ai  orbital 
doubly  occupied. 

Figure  2.  The  nonrigid  intramolecular  rearrangement  for  the  ^E' 
anionic  state  with  the  linear  structures  (^Zg'*')  corresponding  to 
minima  and  the  C2v  structures  (^82)  to  transition  states.  The  reaction 
path  goes  through  Cs  geometries. 


